Scalable uniform construction of highly-conditional quantum gates 
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We present a scalable uniform technique for construction of highly conditional C n -NOT quantum 
gates of trapped ion qubits, such as the Toffoli gate, without using ancilla states and circuits of an 
exorbitant number of concatenated one- and two-qubit gates. Apart from the initial dressing of the 
internal qubit states with vibrational phonons and the final restoration of the phonon ground state, 
our technique requires the application of just a single composite pulse on the target qubit and is 
applicable both in and outside the Lamb-Dicke regime. We design special narrowband composite 
pulses, which suppress all transitions but the conditional transition of the target qubit; moreover, 
these composite pulses significantly improve the spatial addressing selectivity. 

PACS numbers: 03.67.Lx, 03.67.Ac, 32.80.Qk, 42.50.Dv 



I. INTRODUCTION 



One of the most important highly-conditional quan- 
tum gates is the three qubit control-control-NOT gate, 
known as the Toffoli gate, in which the target qubit is 
inverted if both control qubits are in state |1), and is 
left unchanged otherwise. This gate has a central role 
in quantum error correction [l[; moreover, it forms with 
the one-qubit Hadamard gate a universal set of quantum 
gates 0. The more general C"-NOT gates are often used 
in quantum computing, e.g. as oracles in Graver's search 
[H and to simulate quantum walks [IJ. 

The simplest decomposition of the Toffoli gate in the 
circuit model of quantum computation uses six CNOT 
gates [H, or five CNOT gates with an ancilla state [jj. 
Extending this approach to C™-NOT gates with n > 2 
is highly demanding for it requires the ability to con- 
struct, with very high fidelity, many concatenated gates 
of this type. It is hence desirable to seek simpler schemes 
for conditional multi-qubit gates without sequences of 
CNOT gates Recently, Monz et al. have demon- 
strated experimentally the Toffoli gate with 71% fidelity 
with 40 Ca + ions by using a sophisticated sequence of 15 
laser pulses in the Lamb-Dickc (LD) regime; it is not ob- 
vious, however, how this approach can be extended to 
higher C™-NOT gates and outside the LD regime. 

In this paper, we propose a simplified uniform scheme 
for construction of C n -NOT gates of arbitrary order n in 
a linear ion string by using specially designed compos- 
ite pulses. The method does not use ancilla states and 
circuits of concatenated one- and two-qubit gates, and is 
applicable both inside and outside the LD regime. We 
design special narrowband (NB) composite pulses whose 
excitation profiles allow us to drive only transitions be- 
tween a pair of selected collective ionic states, thereby 
manipulating the target qubit in a way, controlled by the 
other qubits. The method allows us to construct C™- 
NOT gates of various orders n with essentially the same 
composite sequences, and hence the same level of com- 
plexity. 
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FIG. 1: Left: Ratios of pulse areas A^+i/A^s with 2 ^ 
v ^ 7, vs the LD parameter r), needed for the Toffoli gate 
when the ion chain is initialized with vq — 5 phonons. Right: 
Excitation profile produced by a three-component BB pulse, 
B3, with composite phases (0 , 0.657T, 0). Note that for a broad 
range of pulse areas around ir the transition probability p for 
B3 is very close to 1, p > 1 — 10~ 4 (horizontal dashed lines). 



II. BASIC STEPS 

Our method begins with the initialization of the string 
of n + 1 ions in the collective vibrational ground state 
\v = 0) 3. We adopt the wavefunction notation |^>)|t>), 
where \ip) = \q1q2 ■ ■ • Qn+i) is the collective internal state 
of the ion qubits, with q^ = or 1, and \v) is the collec- 
tive phonon state. The C™-NOT gate is constructed in 5 
steps. 

Step 1. The ion string is prepared in a common vq- 
phonon Fock state of a selected phonon mode, |0) — > \vo), 
with Vq ^ 71+ 1; this can be achieved with vq alternating 
blue- and red-sideband tt pulses on an ancilla ion. 

Step 2. A three-component broadband (BB) compos- 
ite pulse is applied on each ion on the first blue sideband, 
i.e. 3(n + 1) pulses in total. This sequence is "seen" as 
a 7r pulse by the transitions |0)|i>) <-> + 1), where 
vq — n — 1 ^ v ^ Vq + n, and therefore the respective 
states in these transitions are inverted. This is illustrated 




FIG. 2: Linkage diagram for a system of three ions, in which 
the third ion is addressed with a laser tuned on the second 
blue sideband of a selected vibrational mode. The laser field 
couples states with v and v+2 phonons with coupling strength 
dependent on v and the LD parameter 77, Eq. (JTJ . 



for 3 ions in Fig. [TJ if the transition \0)\v) f> |l)|u + 1) 
with v = 4 phonons "sees" a ir pulse (^4,5 ~ 7r), then 
all transitions with 2 ^ v ^ 7 will "see" (with very 
high fidelity) a ir pulse too. Then the collective inter- 
nal states will be dressed with a different number of 
phonons conditional on the number of qubits in states 
|0) and |1): a collective state with rik qubits in state \k) 
(k = 0, 1) will be dressed with \vo + ri\ — no) phonons. 
For 3 ions states |000)|«o) and 1 1 11) | z^o) are mapped re- 
spectively to states I 111) |i^o + 3) and |000)|do — 3), states 
|100)|uo) and 1 101) |fo) are mapped respectively to states 
|011)ju + 1) and [010)1^0 - 1), etc. Thus the collec- 
tive states |gi 92 ■ ■ • 5n+i)|' y ) group into sets with the same 
number of phonons, and the same total number of inter- 
nal excitations, as shown in Fig. [5J Note that after this 
step the state of each qubit is inverted, \qk) |1 — %)• 
Step 3. A ./V-component composite pulse sequence is 
applied on the target ion; this sequence is the core of our 
method. It must act in such a way that the transition 
|0i0 2 • • • 0„0 n +i) <-> |0i0 2 • • • n l n+ i) "sees" an effective 
7r pulse and is inverted, while all other transitions "see" 
effectively a On or 2n pulse, i.e. remain either unchanged 
(for On) or all acquire the same phase n (for 2n). (It 
turns out that the latter option is easier to realize.) Such 
a discrimination is made possible by step 2, because the 
couplings fl VjV+ 2{t) depend on the LD parameter 77 and 
the number of phonons v, which is different for the dif- 
ferent sets 
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(1) 



where L%(x) is the generalized Lagucrrc polynomial. 
The different couplings produce different pulse areas, 
A v , v +2 = f t . f Qv,v+2(t)dt- Since two adjacent sets differ 
by 2 phonons we have to address the ions on the second 
blue vibrational sideband of a selected mode. 

Consider a string of 3 ions needed for the Toffoli C 2 - 
NOT gate. To perform the Toffoli gate conditional on the 
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FIG. 3: Left: Ratios of pulse areas Ae,s/A2,4 and A^fijA^^ 
vs the LD parameter 77, when the ion chain is initialized with 
Do = 5 phonons. Right: Excitation profile produced by a NB 
pulse N13. The range of gate operation is rj < 0.3, where the 
two pulse area ratios are in the vicinity of 2. 



states of ions 1 and 2 we address only ion 3 on the second 
blue sideband. The laser field then produces four sepa- 
rate two-state systems with different couplings, as shown 
in Fig.[5J We need a composite sequence, which is "seen" 
as a 7r pulse by the transition |000)|uo — 3) -o- |001)|wo — 1) 
but as a 2-7T pulse by the other three transitions. Here the 
dependence of the couplings on the LD parameter 77 is es- 
sential because the variation of 77 allows us to find ranges 
where these conditions arc fulfilled. For an initial phonon 
number vq = 5 the three subsystems in Fig. [5] have cou- 
plings 02,4, ^4,6 and 06,8- The dependence of the pulse 
area ratios ^6,8/^2, 4 and A^^/A^a on 77 are illustrated 
in Fig.[3Jleft). The range of values of 77 where these ratios 
are near 2 are suitable for realization of the Toffoli gate, 
as indicated by the composite-pulse excitation profile in 
Fig. fright). One can vary the initial phonon number 
vq and produce different behavior of the relative pulse 
areas versus 77 according to Eq. ([l}, which allows one to 
shift the 77-range for the Toffoli gate to the experimentally 
most convenient values. 

Step 4- We repeat step 2 [ll|; the phonon number in 
each collective internal state is restored to uo- 

Step 5. Step 1 is applied in a reverse manner; the ion 
chain is restored in its vibrational ground state \v = 0). 

Steps 1-5 produce the transformation 



|V>|0> -> [C' l -NOT|V>]|0), 



(2) 



with 8(71+1) single pulses and a composite sequence. The 
number of pulses in the composite sequence from step 3 
can vary, e.g., from 5 for C 2 -NOT (Toffoli) gate to 13 for 
C 6 -NOT gate, as shown below. Because the number of 
pulses in the composite sequence is about 2n, the total 
number of pulses is about lOn. 

Our method generalizes the original idea of Monroe 
et al. [ijj who constructed a CNOT gate by exploiting 
the D-dependence of the LD parameter 77 to select such 
a "magic" value of 77, for which the ratio of the two rel- 
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standard NB pulses 
N 5 (1.160; 0.580) 
N 9 (1.130; 0.820; 0.110; 1.390) 
Ni 3 (1.270; 0.440; 1.020; 0.770; 1.850; 1.730) 
N17 (1.600; 0.550; 1.090; 0.890; 0.620; 1.540; 0.150; 1.570) 
N 2 i (1.070; 0.920; 0.130; 1.830; 1.160; 0.720; 0.100; 1.520; 

0.810; 1.950) 

N 25 (1.750; 0.380; 1.420; 0.710; 1.070; 0.910; 0.780; 1.470; 

0.550, 1.740, 0.160, 1.650) 

optimized NB pulses 
Ng (1.190; 0.630) 
N§ (1.157; 0.888; 0.218; 1.529) 

N 13 (0.585; 1.352; 0.914; 1.186; 0.020; 0.089) 



TABLE I: Phases (<j>2', (j>i\ <^4> ■ • ■ i 0m+l) (in units ir) for stan- 
dard (Njv) and optimized (N%) NB sequences of N = 2m + 1 
pulses of area A: AoA^A^ ■ ■ ■ A (t>m+1 ■ ■ ■ A^A^Aq. 

evant couplings equals a rational odd/even or even/odd 
number. This approach cannot be extended with a single 
pulse to higher gates as they require specific values of the 
ratios of three or more couplings. We can satisfy these 
conditions in various 77-ranges because composite pulses 
can modify the excitation profile in any desired manner. 



III. COMPOSITE SEQUENCES FOR C"-NOT 
GATES 

The technique of composite pulses was introduced in 
nuclear magnetic resonance (NMR) [l3|, [l4[ as a powerful 
tool for control of spins by magnetic fields. A compos- 
ite pulse is a train of TV pulses with well-defined relative 
phases <j>k (k = 1,2,..., TV), which are used as control 
parameters in order to compensate the imperfections of 
a single pulse and/or to shape up the excitation pro- 
file in a desired manner. We have designed special com- 
posite pulses for the C"-NOT gates by using a recently 
developed simple method [HI, [l6|. We use equal pulse 
areas A and an odd number of pulses, N = 2m + 1, al- 
though these restrictions are not essential. We consider 
symmetric "anagram" composite sequences, with phases 
4>k = 4>N+i-k (k = 1,2, ...,m); this condition pro- 
duces symmetric excitation profiles. Because the over- 
all phase of the composite sequence is irrelevant we set 
4>i = 4>n — 0; hence we are left with m different phases. 

Three families of composite pulse sequences are par- 
ticularly important: broadband (BB), narrowband (NB) 
and passband (PB) [TJ]. The BB pulses stabilize pop- 
ulation inversion to values p « 1 around the pulse area 
7r (flat-top excitation profile) . The NB pulses stabilize p 
to values p m around the pulse area Ott (or 27r) (flat 
bottom). The PB pulses stabilize p both to values p « 1 
around area n and to values p « around area On (or 
27r) (flat-top and flat-bottom). 



The composite phases for the BB sequence used in 
steps 2 and 4 are (0, 0.657T, 0). Other examples of BB 
sequences can be found in 1 1 €31 ] . 

In step 3 all but one of the couplings must fall in the 
region where p « 0, hence stabilization is needed around 
this value; the remaining coupling must fall in the range 
where p « 1, i.e. near it. This condition suggests to use 
NB pulses because of their flat-bottom excitation profiles 
around areas Oir and 2ir. None of the existing NB com- 
posite pulses, however, produces excitation profiles with 
sufficiently wide bottoms to satisfy these conditions with 
high fidelity; for this reason we construct here new NB 
composite sequences. We derive the NB phases from the 
conditions 

[d k A u[^] A ^ = (k = 2,4, ...,2m), (3) 

with d A = d k /dA k ; the skipped derivatives vanish identi- 
cally for "anagram" sequences. \j( N > is the full propaga- 
tor, UW = U O U 02 U 03 • • • U 0m+1 • • • U 03 U 02 U O , with 

u _ cos(A/2) ie-^sin(A/2) (4) 
~ ie^sm(A/2) cos(A/2) [ ' 

Examples of NB sequences suitable for C™-NOT gates are 
presented in Table HI A typical NB excitation profile is 
plotted in the right frame of Fig. [3] These NB sequences 
are "all-purpose" ones, i.e. they are suited whenever an 
increased selectivity of excitation, or suppression of un- 
wanted excitation is needed, e. g. f or local addressing in a 
lattice of closely spaced qubits [15| . We have constructed 
special NB sequences which optimize the performance of 
the C"-NOT gates; they are listed in Table U as well. 
Their phases are obtained by starting from a particu- 
lar "all-purpose" NB pulse and minimizing the error in 
the high-fidelity region of the respective C n -NOT gate in 
the (77, A) control landscape. We can always accommo- 
date more couplings (that have to be suppressed) within 
the no-transition bottom of the excitation profile around 
27T (and 07r) by adding more pulses to the composite se- 
quence in order to broaden this bottom. 

After the NB sequence the state of the target ion is 
changed as |0) — » c~^|l) and |1) -4 c^|0), where 

m 

^ = (-1)'" (tt/2 + </wi) - 2^(-l) fe ^- (5) 

k=2 

This phase can be compensated, if necessary, with an 
additional Stark pulse focused on the target ion. 

IV. SIMULATION OF C' -NOT GATES 

We have simulated numerically different C n -NOT 
gates C n by using our approach presented above. The 
fidelity of the simulated gate S„ is defined as the uni- 
form average over infinitely many random states 

F = mean^j | (■(/;|Sj l C„|V') | 2 - Figure 2] shows the fidelity of 
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FIG. 4: Simulated fidelity of the Toffoli gate with NB compos- 
ite sequences from Table [I] for initial phonon number vo — 4 
(solid lines) and 5 (dashed lines) vs. the LD parameter rj and 
the pulse area A. Left: standard NB sequences N5 (top) and 
Ng (bottom). Right: optimized NB sequences N5 (top) and 
N9 (bottom). The area is divided by A m i n = A V0 -3 yV0 -i, that 
is the (smallest) pulse area which has to be "seen" as a it pulse 
by the transition |000>|u -3) o |001)|«o - 1), cf. Figs, [land 
1 



the TofToli gate constructed with standard (left) and opti- 
mized (right) NB sequences from Table [I] for 5 (top) and 
9 (bottom) ingredient pulses. Fidelity above 99%, and 
even 99.9%, can be obtained with just 5 pulses. Higher 
fidelity can be obtained with longer sequences, e.g. over 
99.99% with the Ng pulse. The high-fidelity 77 ranges can 
be varied by varying the initial phonon number Vq, as 
evident from the examples with vq = 4 and 5 in Fig. [4] 
a larger vq shifts the high-fidelity 77-range toward 77 = 0. 

We have simulated several higher-order C™-NOT gates 
with n = 3 to 6 using the NB sequences from Table [H 
the fidelities for the optimized NB pulse N° 3 are shown in 
Fig. [5J High fidelity can be achieved in various ranges of 
77 by changing the initial phonon number vq because the 
ratios of the couplings depend on both 77 and vq; larger 
Vq push the high-fidelity range again toward smaller 77. 
In this manner, one can adjust the operation of the C n - 
NOT gate to a range of values of 77, which is most suit- 
able experimentally. Fidelity can always be increased, 
and the high-fidelity ranges can be expanded, by adding 
more pulses to the composite sequence, which will allow 
enhanced optimization. It is remarkable that the same 
composite pulse NJ 3 implements all gates up to n = 6 
efficiently, in similar ranges of 77 and A. We have checked 
that this same pulse N° 3 can be used with over 99% fi- 
delity for C"-NOT gates up to n = 10. The lower gates 
with n = 3 and 4 can be implemented efficiently also with 
the shorter pulse Ng (not shown), which has been used 
in Fig. [3J This uniformity is strikingly different from the 
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FIG. 5: Simulated fidelity of C"-NOT gates with n = 3 
to 6 with the optimized NB sequence N° 3 from Table [J vs. 
the LD parameter 77 and the pulse area A. The area is di- 
vided by Ami n — Auo-n-iyoo— n+ij that is the (smallest) pulse 
area which has to be "seen" as a 1 pulse by the transition 
|00---00)|i>o-n-l) |00---01» -n + l). Note that the 
same composite sequence NJ 3 implements all gates efficiently. 



circuit model, in which the complexity of implementation 
increases rapidly with n, which makes the demonstration 
of most algorithms still impossible with current technol- 
ogy 0- We have verified that, for up to n = 20, NB 
composite sequences of about 2n ingredient pulses are 
sufficient for efficient construction of the C n -NOT gate. 



V. DISCUSSION AND CONCLUSION 

The proposed method offers a conceptually simple and 
scalable implementation of C"-NOT gates of trapped- 
ion qubits of arbitrary order n, as it requires, besides the 
usual initial dressing and the final undressing of the col- 
lective qubit states with phonons, just a single NB pulse 
applied on the target ion qubit. This method uses the 
basic physical notion of destructive interference of un- 
wanted transitions rather than quantum circuits of a vast 
number of concatenated one- and two-qubit gates. The 
ensuing simplicity and universality make our method eas- 
ily scalable to arbitrary C"-NOT gates with essentially 
the same level of complexity. The constructed NB se- 
quences produce C™-NOT gates with very high fidelity 
using relatively few ingredient pulses. We have found 
that a NB composite sequence of about 2n pulses suffices 
for a C™-NOT gate; including the dressing and undress- 
ing steps, the total number of pulses is about lOn. This 
makes possible the creation of C™-NOT gates for up to 
n = 15 — 20 ion qubits with the existing ion trap technol- 
ogy [13 ■ We point out that the NB composite sequences 
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are very convenient in another respect: they eliminate 
spatial imperfections of local addressing (l5j . which are 
often a limiting factor in experiments |7|. 

We note that the use of blue-sideband pulses is not 
mandatory; one can use red-sideband pulses too. One can 
also replace the NB composite sequences, which minimize 
the number of ingredient pulses, by PB sequences, which 
require a few more ingredient pulses but provide in return 
greater robustness against variations in the pulse area. 

We point out in conclusion that the initial dressing of 
the qubit states with phonons and the final restoration of 



the phonon ground state, each of which requires about 4n 
pulses, can be optimized for each C™-NOT gate to fewer 
pulses, about n for each of the dressing and undressing 
steps [18(. The procedure described above is, however, 
universal (applicable to arbitrary n) and flexible, for it 
applies to large ranges of values of the LD parameter ?y; 
the latter feature allows one to operate outside the LD 
regime and hence speed up the gate operation. 
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sion project FASTQUAST and the Bulgarian NSF grants 
D002-90/08 and DMU02-19/09. 
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